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Abstract. In this paper, ten CME events viewed by the 
STEREO twin spacecraft are analyzed to study the deflections 
of CMEs during their propagation in the corona. Based on the 
three-dimensional information of the CMEs derived by the gradu- 
ated cylindrical shell (CCS) model | Thernisien it is 
found that the propagation directions of eight CMEs had changed. 
By applying the theoretical method proposed hy \Shen et al.\ |2011| 
to all the CMEs, we found that the deflections are consistent, in 
strength and direction, with the gradient of the magnetic energy 
density. There is a positive correlation between the deflection 
rate and the strength of the magnetic energy density gradient 
and a weak anti-correlation between the deflection rate and the 
CME speed. Our results suggest that the deflections of CMEs are 
mainly controlled by the background magnetic fleld and can be 
quantitatively described by the magnetic energy density gradient 
(MEDG) model. 

1 Introduction 

Corona mass ejections (CMEs) are large scale eruptions 
from the solar surface and act as one of the primary drivers 
of space weather phenomena, such as geomagnetic storms, 
solar energetic particle events, etc. The deflections of CMEs, 



which were first reported bv \MacOueen et al\ jl986t in the 
Skylab epoch (1973-1974), are one of the factors influenc- 
ing the geoeffectiveness of CMEs. A stat istical study about 
CME deflections was made by Crem ades and BothmeA 
|2004l ]- They identified the source regions of 124 structured 
CME events observed by the Large Angle and Spectrometric 
Coronagraph (LASCO) on board the Solar and Heliospheric 
Observatory (SOHO) during 1996-2002 and compared the 
source re gions of CMEs with the ir central position angles 
(CPAs). \Cremades and BothmeA j2004 found that there 



was a systematic deflection by 20 degrees to lower latitudes 
only in activity-minimum years (1996-1998) and no system- 
atic trend nor deflection during the years 1999-2 002. The re- 
sult w as further confirmed by the recent work o AWano et al\ 
[2011^ . in w hich all the LASCO CM Es during 1997-1998 were 
examined. \Cremades et al\ |2006[ | found a good correspon- 
dence between the deflection of CMEs and the total area of 
coronal holes (CHs). They suggested that the neighboring 
CHs affect the outward evolution of CMEs near the Sun and 
cause such deflections. \Sh en et al\ |201l[ | analyzed the de- 
flection of the 8 October 2007 CME in the meridian plane in 
much more details. They showed strong evidence that the 
trajectory of the CME was influenced by the background 
magnetic field, and the CME tends to defiect to the region 
with lower magnetic energy density. 

Note that all the CME defiections studied above are in the 
latitudinal direction. The CME defiections on a spherical 
surface, i.e., in both latitudinal and longitudinal directions, 
still remain unclear due to the presence of projection effect. 
Even so, the CME defiect ion in longitude was suggested by 
some researchers. For example, the longitudinal defiections 
of CMEs as they propagate in the interpla netary space were 
first proposed hw iwano etal\ j2004l . l2006at . Such defiections 
can explain the asymmetrical east-west distri bution of the 
sourc e locations of the geoeflFectiye halo CME s [Wano et all 
12002^ . ^Govalsw amv et al\ j2004l . l2005l . 12009^ also suggested 
that CMEs could be defiected away from the Sun-Earth line 
by the associated coronal holes. They use such defiections to 
explain the existence of the 'driver less' shocks, which were 
observed near the Earth but without their drivers, the inter- 
planetary coronal mass ejections (ICME). 
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Since the successful launch of the Sol ar TErrestria l REla- 
tions Observatory (STEREO) mission [Kaiser et all . [2008], 
the three-dimensional (3-D) information of CMEs is more 
or less revealed in observa tions with the aid of various 
reconstruction models [e.^., Thernisien et all l200d l2009l : 
\Luaaz et al\\201(i : lLiu et a/.Ll2QlQal lbt. STEREO consists of 



two identical satellites. It has provided for the first time the 
observations of the Sun from dual vantage points. Based 
on the STEREO observations, some CME events with an 
obvious defl ection in the latitud i nal direction have been re- 
ported [e.g., \Kilvua et all . l2009l : \Shen et ad 1201 ll] . and the 
possible deflections o f CMEs in the ecliptic plane ha v e also 
been discussed [ e.^.. \ Liu et all l2010bl : \Luaaz et all |201Q| : 



\Poomvises et all 1201' 

In this paper, we have comprehensively studied the CME 
deflections in the corona in both latitudinal and longitudi- 
nal directions for ten CME events viewed by the STEREO 
twin spacecraft. The data and the method we used will be 
introduced in the next section. In Section [3l four cases will 
be selected to show different types of deflections, in which 
their 3-D trajectories and a comparison between the deflec- 
tions and the magnetic energy density distributions will be 
presented. In Section [H statistical studies on the deflection 
and its correlation with the magnetic energy density will be 
presented. Finally, we will give the conclusions and make 
some discussion in Section O 

2 Data and Method 



CCS. 

A set of reasonable initial values of the parameters is help- 
ful to get the best fitting of the CME images. Observations 
of the CME source region on the solar surface were used to 
constrain the initial values of longitude 0, latitude and tilt 
angle 7 if any. The tilt angle can be estimated according 
to the CME-associated filament (or the polarity inversion 
line, PIL, if no filament observed) because it is believed that 
a CME is a flux rope surrounding its associated filament 
and standing above the PIL. The rest of the parameters are 
set by comparing the CCS flux rope to the CME shape ob- 
served simultaneously by both STA and STB. In practice, 
we find that the tilt angle and half angle are insensitive to 
the fitting results. Therefore we fix them to a certain rea- 
sonable value for the whole CME evolution process by trial 
and error. It should be noted that fixing the tilt angle in- 
dicates a CME without rotation, which may not be true for 



many CME events 


Lynch et all l2009l: \Mdstl et a/.. 2008: 


IShiota et al, 2010l: 


Torok and Klieni, l2003l: \Wana et al.. 


2006bl; Yurchvshvn et a/.,l2007ll. However, we find that the 



change of the tilt angle will not significantly affect the de- 
rived directions of CMEs as long as the CCS flux rope 
matches the observed CME shape in both the STA and STB 
images (see the discussion in Section [5|). Thus, in this study 
a fixed tilt angle is acceptable. Besides, not all the CMEs 
in our sample have available observations of their source re- 
gions. For such events, we just compare the CCS flux rope 
with the observed CME shape to estimate the parameters. 



2.1 Three-dimensional information of CMEs 2.2 Coronal magnetic field 



The observations from the CORl and the C0R2 instru- 
ments of the Sun Earth Conne ction Coronal and Helio- 
spheric Investigation (SECCHI) [ Howard et all |2008[ | suite 
on board the STEREO A and B spacecraft are used to learn 
about the evolutions of CMEs in the corona. The CORl in- 
struments observe the corona from 1.4-4.0 Rs and the C0R2 
instruments observe the corona from 2.5-15.0 Rs- In this 
paper, these observations were used to obtain the 3-D infor- 
mation of CME during its propagation in the corona. The 
observations from the SECCHI/EUVI a nd the Michelson 
Dopp ler Imager (MDI) on board SOHO \Brueckner et all 
Il995[ | are used to identify the CMEs' source regions on the 
solar surface. 

To obtain the 3-D geometry and therefore the trajectory 
of a CME, the graduated cyhndrical sh ell (CCS) model de- 
veloped hY \Thernisien et al\ [20061 . 12009^ was applied to both 
the projected two-dimensional images from the STEREO-A 
(STA) and the STEREO-B (STB) spacecraft. In that model, 
CMEs are assumed to have a flux rope-like structure. The 
GCS model has nine f ree parameters (refer to Table 1 of 
\Thernisien et al\ [20061 ]). Six of them determine the CME's 
shape projected on the plane of the sky. These parameters, 
referred to as geometric parameters, are the longitude '0' 
and latitude height '/i' (the height of the legs, or '/i/', 
the height of the leading edge), aspect ratio 'a^', tilt angle 
'7' with respect to the equator, and half angular width 'a' 
between the flux rope legs. The other three parameters, 
specifying the electron density distribution at the shell, are 
the electron density factor W^', Gaussian width ^ a trailing^ of 
the density profile in the interior of the GCS and Gaussian 
width ^aieading^ of the density profile at the exterior of the 



It is believed that the magnetic energy is dominant in 
the corona. Previous studies have suggested that the 
CME deflection can be qualitatively interpreted as the con- 
straint of the ambient ma gnetic structure, e.g., coronal holes 
[Govalswamv et oil . |2004| ]. Our recent study of the 8 Oc- 
tober 2007 CME showed that the behavior of the CME's 
latitudinal deflection can be quantitatively described by a 
theoretical method, in which the direction and magnitude of 
the deflection are well consistent with the gradient (with the 
conventional minus sign in front) of magnetic energy density, 
< — V(^;^) >, where the angle brackets mean the average 
over the region occupied by the CME \Shen et 'all , l201l[ | . In 
that work, however, the deflection in the latitudinal direc- 
tion of only one CME was studied. Thus, we will test the 
method with more CME events to check if it is also applica- 
ble to other CMEs and to the deflections in other directions. 

In the method, the 3-D magnetic field of the 
corona is key information, and is extrapolated from the 
SOHO/MDI photospheric magnetic synoptic charts by the 
current-sheet source - surfac e (CSSS) model developed by 
\Zhao and Hoeksemol [l995l] . The CSSS model is a devel- 
opment of the potential-field sourc e- surface (PFSS) niodel , 
and was used in our pre vious work [Shen etal. 1, 12OO7I . I2OIII : 
I Wang and ZhanX l2007l ]. The magnetic synoptic chart is 
created from the MDI daily magnetograms over a quasi-27- 
day solar rotation. It cannot reflect the state of the photo- 
sphere right at the time of a CME taking place. However, 
what we are interested in is the large scale coronal magnetic 
field, which probably ch anges little during a solar rotation 
\Ness and VF^z/coxl . Il964l ]. Under this consideration, the syn- 
optic chart may be treated as a good approximation to the 
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(a) 



(b) 



Figure 1: The source region of the 12 December 2008 CME observed by the STEREO/EUVI and the SOHO/MDI. Panel 
(a) shows the STEREO/EUVI 171 A image superposed by the contours of the SOHO/MDI magnetogram taken two days 
before the event when the source region (marked by the square box) was visible to the SOHO. The STEREO/EUVI image 
was rotated to match the angle of view of the SOHO/MDI. Panel (b) shows the zoomed-in image of the source region. 




A:2008-1 2-1 2107:05:00.007 



A:2008-1 2-1 2107:05:00.007 



A:2008-1 2-1 2107:05:00.007 



3 



B:2008-1 2-1 2107:05:35.865 



B:2008-1 2-1 2107:05:35.865 



B:2008-1 2-1 2107:05:35.865 



Figure 2: The fitting example of CME-1. From (a) to (c): The original CME images, the modeled wireframe images which 
overlay on the CME images, and the relative brightness derived from the CCS model. The top and bottom panels present 
the results based on the STA and STB data, respectively. The arrows with numbers indicate the common points of the 
CME between the original images and the derived brightness images. 
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real photospheric magnet ogr am over the full solar surface. 
In this study, the magnetic synoptic charts with spatial res- 
olution of 360° X 180° are used. To get the best extrapolation 
results, the order of the harmonic coefficients is chosen to be 
125. Once the coronal magnetic field is extrapolated, the av- 
erage gradient of the magnetic energy density can be easily 
calculated for any CMEs of interest. 

3 Observations and model analyses of four 
cases 

Before we show the statistical results of ten CMEs, in this 
section four different types of CMEs are selected to investi- 
gate in detail the deflection behaviors and their relationship 
with the gradient of the coronal magnetic energy density. 
The first case is the 12 December 2008 event (CME-1), which 
deflected in the latitudinal direction. The second case is the 
9 April 2008 event (CME-2), which deflected in the longi- 
tudinal direction. The third case is a CME erupting on 16 
November 2007 (CME-3), in which a deflection in both lat- 
itudinal and longitudinal directions was obvious. The last 
case is the 3 November 2008 event (CME-4), which did not 
show an evident deflection. 

3.1 The 12 December 2008 Event (CME-1) 

This CME first appeared in the field of view (FOV) of 
the STA/CORl and the STB/CORl at about 05:35 UT on 
12 December 2008. To fit the CME with the GCS model, 
it is required that the CME almost fully appeared in the 
FOV of the coronagraph. Thus the first and the last CORl 
image pairs we selected are taken at 05:35 UT and 07:35 UT, 
respectively, during which there are 13 image pairs (or data 
points) with a cadence of 10 min. Similarly, in C0R2 FOV, 
the first and last image pairs are taken at 09:52 UT and 
14:52 UT, respectively, and there are 11 image pairs with a 
cadence of 30 min. 

This CME was associated with a filament, which erupted 
at about 04:00 UT on 12 December 2008. The square area 
in Figure [ija) represents the source region of the CME by 
combining the STEREO_A/EUVI and SOHO/MDI observa- 
tions. The MDI magnetogram is superposed on the EUVI 
171 A image as indicated by the red (positive) and green 
(negative) contours, and the EUVI image is rotated to match 
the time and the vantage point of the MDI data. The black 
curve in the square displays the filament. This filament ex- 
tended over a long and narrow region, from about 77° to 
106° in longitude and about 29° to 49° in latitude under the 
Carrington coordinate system. 

Figure mb) shows the zoomed-in image of the CME source 
region. As has been mentioned in the last section, the tilt 
angle and the half angle are fixed to a certain value for the 
whole evolution process of the CME. After applying a trial 
and error method, we find that, with the value of —15° for 
tilt angle and 14° for half angle, the GCS model can reach 
the best fitting result by a visual judgement. 

Further, we fit the observed CME shape with the GCS 
model for each image pair. Figure [2] shows the sample of 
the fitting result of the CME recorded at 07:05 UT. Figures 
[2jbl) and[2l^b2) present the wireframe of the GCS flux rope 
which is overlaid on the original images. From these two 



Table 1: The fitted free parameters of 
CME derived by the GCS model with 
and the half angle of 14°. 



the 12 December 2008 
the tilt angle of —15° 
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images, we found that the shapes of the CME are both con- 
sistent with the wireframe. Figures [2]^cl) and[2jc2) present 
the relative brightness images derived from the GCS model. 
By comparing Figures [2jc) with[2l^a), we find that they are 
quite similar. The arrows marked in these panels denote 
some common points: the bright features in Figure [2j a) are 
also bright in Figure [2jc) (arrows 1 and 3), and the darker 
features in Figure [2ja) also look darker in Figure [2jc) (ar- 
rows 2 and 4). These results indicate that the GCS model 
not only fits the projected shape of the CME well, but also 
could explain the relative brightness of the CME. Based on 
the above analysis, we are quite confident that the derived 
parameters should well reflect the 3-D geometry of the CME. 

After all the 24 image pairs are processed, the CME tra- 
jectory is obtained. All the fitted free parameters are listed 
in Table [T] The first column gives the time when the CME 
was recorded by the STEREO/SECCHI instrument. The 
next two columns give the longitude under the Carrington 
coordinate system (^0c ' ) and the Stonyhurst coordinate sys- 
tem ('(/)s') I Thomvsoj\ . I2OO6I ] , respectively. The next three 
columns give the other geometric parameters: latitude 
height '/i/', ratio 

Figure [31^ a) shows the height-time plot of the CME. Both 
linear and quadratic fittings are applied to the measure- 
ments. It is found that this CME was propagating outward 
with a speed of about 275.6 km and an acceleration of 
about 12 m . Figure [3jb) and [3jc) present the Stony- 
hurst longitude and latitude as a function of the height, re- 
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Figure 3: The kinetic evolution of the CME-1 event. Panels 
(a) to (d) show the height-time, Stonyhurst longitude-height, 
latitude- height, and deflection rate- height plots, respec- 
tively. The error bars in the first three pa nels are got from 
the 10 % decrease of the merit function (see \Thernisien et all 
[2009^ for details). 



spectively. The Stonyhurst longitude changed around the 
value of 5° with ^ 4° variation. Considering the error in 
our fitting process, this CME did not manifest an evident 
deflection in the longitude. But its latitude shows a clear 
variation from about 30° to 10°, which suggests that the 
CME experienced an evident deflection from high latitude 
to low latitude. 

Further, we define the deflection rate as Aa/A/i, where 
Aa is the deflection angle (both latitudinal and longitudi- 
nal deflection are taken into account). For events with more 
than 10 data points, the deflection rate at any data point 
is calculated by fitting the longitude and the latitude with 
height of neighboring five data points. For events with less 
than 10 data points, a fitting procedure over three neigh- 
boring data points is used. The variation of the deflection 
rate with the height is shown in Figure [Sj^d). It is clear that 
the deflection rate decreases quickly as the height increases. 
The main deflection of the CME occurred in the range below 
about 8 Rs. When the CME's leading edge exceeded 8 Rs, 
the deflection became insignificant. This event has been pre- 



viously studied by some other re searchers \e.Q:.\Bvrne et al 
2010; Davis et al, 2009| ; [lIu et a /.. 2010h: IPoomvises et al 



2010l : iLuaaz et a/.ll2QTQ| . The results we obtained here are 
consistent with their results. 

On the other hand, the magnetic field energy density dis- 
tribution at corresponding altitude for every data points is 



calculated. The SOHO/MDI magnetic synoptic chart of 
the 2077 Carrington rotation which begins at 07:00 UT, 20 
November 2008 and ends at 14:38 UT, 17 December 2008 
is used as the bottom boundary for the CSSS model. Fig- 
ured! shows the distributions of the magnetic energy density 
at different altitudes. The red curves indicate the position 
of the heliospheric current sheet (HCS, only marked in the 
panels with the altitude larger than 2.6 Rs, where the coro- 
nal magnetic field is open). The yellow asterisk marks the 
projected location of the CME leading edge on the Carring- 
ton map, and the cyan ellipse indicates the boundary of the 
CME in projection. The average value of the gradient of the 
magnetic energy density in the ellipse is marked by the red 
arrow, and the corresponding CME deflection is marked by 
the green arrow. The length of the arrows indicates the rel- 
ative strength of the gradient and the deflection rate. The 
length is scaled by comparing its strength with all the data 
points of the ten CMEs. From the figure, it can be seen that 
the CME deflection is consistent well with the gradient of 
the magnetic energy density in both strength and direction, 
which roughly points from high latitude to low latitude. As 
a consequence, the CME leading edge was getting closer to 
the HCS during its propagation. 

3.2 The 9 April 2008 Event (CME-2) 

This CME first appeared in the FOVs of the STA/CORl 
and STB/CORl at about 10:45 UT on 9 April 2008. To 
guarantee that the CME almost fully appeared in the FOV, 
the first and last image pairs of CORl data were taken at 
10:45 UT and 11:25 UT, respectively. There are five image 
pairs during the interval. The first and last images of C0R2 
data are taken at 13:22 UT and 14:52 UT, respectively, and 
a total of four image pairs are selected. 

The CME was associated with an eruptive filament which 
erupted at about 09:21 UT on 9 April 2008 seen from the 
STB. Figure [5]^ a) represents the combined image of the 
STEREO_A/EUVI 171 A data and the SOHO/MDI data. 
The square box denotes the source region of the CME. The 
filament marked by the black line located from about 198° 
to 206° in Carrington longitude and about —13° to —15° in 
latitude. Figure [5jb) shows the detailed image of the source 
region. 



Table 2: The fitted free parameters which derived by the 
model with the tilt angle of 8° and the half angle of 11° of 
the 9 April 2008 CME. 
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Figure 4: The comparison between the gradient of the magnetic energy density and the deflection of CME-1. The 
magnetic field energy density in gray scale in each panel is calculated based on the extrapolated coronal magnetic field at 
the corresponding altitude. The unit of the color bar of is J km~^ in logarithm. The projected leading edge of the CME is 
indicated on the Carrington map by the yellow asterisk. The deflection and the gradient are represented by the green and 
red arrows, respectively. The lengths of the green and red arrows indicate the deflection rate and the relative strength of 
the gradient, respectively. The red curves indicate the heliospheric current sheet, which appears above about 2.6 Rs where 
all the coronal magnetic field lines open. 
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MDI 3-A>1^2008 01:34:27.00 




(a) (b) 

Figure 5: The source region of the 9 April 2008 CME observed by the STEREO/EUVI and the SOHO/MDI. Panel (a) 
shows the STEREO/EUVI 171 A image superposed by the contours of the SOHO/MDI magnetogram taken two days 
before the event when the source region (marked by the square box) was visible to the SOHO. The STEREO/EUVI image 
was rotated to match the angle of view of the SOHO/MDI. Panel (b) shows the zoomed-in image of the source region. 



(a1) 




A:2008-04-09T1 3:52:00.004 



(a2) 



(b1) 



A:2008-04-09T1 3:52:00.004 



(b2) 



(c1) 




A:2008-04-09T1 3:52:00.004 



(c2) 




Figure 6: The fitting example of CME-2. From (a) to (c): The original CME images, the modeled wireframe images which 
overlays on the CME images, and the relative brightness derived from the CCS model. The top and bottom panels present 
the results based on the STA and STB data, respectively. The arrows with numbers indicate the common points of the 
CME between the original images and the derived brightness images. 
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Similar to the fitting procedure applied to CME-1, we fix 
the tilt angle and half angle to 8° and 11°, respectively, by 
trial and error. Then we fit the CME shapes for each image 
pair with the GCS model. Figure [6] shows the sample of the 
fitting result of the CME which was recorded at 13:52 UT. 
The wireframe of the model matches well with the CME 
shapes viewed in both STA and STB spacecraft (Figures 
[6l^a) and[6l^b)). Figures [Gj^cl) and [6l^c2) present the rela- 
tive brightness of the CME derived from the GCS model. 
They are quite similar with the observed bright structure 
in Figures ini^al) and[6l^a2). The arrows in the figure mark 
some example common points between the modeled relative 
brightness images and the observed images. Thus, we be- 
lieve that the 3-D geometry of the CME is reproduced by 
the GCS model. Table [2] lists the other parameters derived 
by the model with the tilt angle of 8° and the half angle of 
11° of all the 9 image pairs of the CME. 
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Figure 7: The kinetic evolution of the CME-2 event. Panels 
(a) to (d) show the height-time, longitude-height, latitude- 
height, and deflection rate- height curves, respectively. The 
error bars in the first three panels ar e obtained from the 10% 
decrease of the merit function (see \ Thernisien et all [2009l | 
for details). 

Figure [TJa) shows the height-time plot of the CME. The 
average speed of the CME is 476.3 km s~^, and the average 
acceleration is 20 m s~^. The variations of the Stonyhurst 
longitude and latitude of the CME are shown in Figures 
\7lh) and [He). Different from CME-1, this CME manifested 
a weak deflection in the longitudinal direction, but no ob- 
vious deflection in the latitudinal direction. Its longitude 
systematically changed by about 16° from ^ 97° to ~ 113° 
though the errors are large. The deflection rate of the CME 



is presented in Figure [Tijd). For CME-2, the fitting of the 
longitude and latitude with height over three neighboring 
data points is used to calculated the deflection rate as there 
are just nine data points in total. Similar to CME-1, the 
deflection mainly occurred at the low altitude where the de- 
flection rate is as large as ~ 10° /Rs, and it quickly decreased 
to about 1° /Rs beyond ?^3-5 Rs. 

The deflections of all the data points and the magnetic 
field energy density distributions of CME-2 are compared 
as shown in Figure [HI The synoptic chart of the 2068 Car- 
rington rotation, which begins at 01:18 UT, 20 March 2008 
and ends at 08:09 UT, 16 April 2008 is used as the bottom 
boundary of the CSSS model. From the figure, it can be seen 
that the gradient direction of the magnetic energy density is 
mainly aligned to the deflection at a low altitude. While at 
a higher altitude, the angle between the two arrows becomes 
bigger. However, the gradient of the magnetic energy den- 
sity and the deflection rate both decreased to quite low levels 
at the high altitude. Same as CME-1, the CME leading edge 
also approached close to the HCS during the propagation. 

3.3 The 16 November 2007 Event (CME-3) 

CME-3 first appeared in the CORl FOV at about 07:35 
UT on 16 November 2007. The first and last image pairs of 
the CORl data were selected at 09:35 UT and 10:35 UT, re- 
spectively. During the interval, there are seven data points. 
The first and the last images of the C0R2 data are taken 
at 13:52 UT and 15:22 UT, respectively, and a total of four 
data points are selected. 

Although CME-3 erupted from the front side of the solar 
surface, no significant surface activity was found. Therefore, 
all the parameters were obtained by the image fitting. For 
this event, the optimized tilt angle and half angle are —25° 
and 8°, respectively, when the GCS model can reach the 
best fitting of the CME images. As shown in Figure [T2l the 
GCS flux rope can fit the CME fairly well. Table [3] lists the 
parameters of the event. 



Table 3: The fitted free parameters which derived by the 
model with the tilt angle of —25° and the half angle of 8° of 
the 16 November 2007 CME. 



Time (^c 
[UT] [deg] 



[deg] [deg] 



hf/Rs 



CORl 



09 


35 


303 


1 


101.1 


-22 


2 


3.04 





25 


09 


45 


306 


9 


105.0 


-22 


1 


3.21 





25 


09 


55 


310 


2 


108.3 


-21 


6 


3.29 





25 


10 


05 


310 


7 


108.9 


-21 


2 


3.36 





26 


10 


15 


311 


6 


110.0 


-21 


2 


3.43 





26 


10 


25 


312 


8 


111.3 


-21 





3.56 





26 


10 


35 


314 


2 


112.7 


-20 





3.71 





26 










C0R2 










13 


52 


318 


1 


118.5 


-14 


3 


8.71 





26 


14 


22 


321 


8 


122.4 


-14 


7 


9.73 





27 


14 


52 


322 


9 


123.8 


-14 


2 


10.93 





27 


15 


22 


322 


4 


123.6 


-13 


6 


11.97 





27 



Figure [TOi; a) shows the height-time plot of the CME. The 
average speed of the CME is 255.0 km s~^, and the average 
acceleration is 19 m s~^. The variations of the Stonyhurst 



8 



no I 

0^1 



20080409-2.52R, 



I 

I 0,0 



051 

uok: 





4.42 




4.34 




4.25 




4.12 






1 









20080409-2,78R, 






4.14 




4.06 




3.96 




3.84 




3.66 


1 


3.36 




1.21 



20080409_3,0JR, 




20080409_5,77R, 




20080409-6,57R, 




20080409-9.98R, 




Figure 8: Same as the Figure H but for the 9 April 2008 CME. 





A:2007-1 1 -1 6109:55:00.006 A:2007-1 1 -1 6109:55:00.006. 
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Figure 9: The fitting example of CME-3. From (a) to (c): The original CME images, the modeled wireframe images which 
overlays on the CME images, and the relative brightness derived from the CCS model. The top and bottom panels present 
the results based on the STA and STB data, respectively. 
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Figure 10: The kinetic evolution of the 16 November 2007 
event. 



longitude and latitude of the CME are shown in Figures 
llOf b) and llOf c). Different from the above two events, the 
CME manifested an evident deflection in both longitude and 
latitude. The Stony hurst longitude systematically changed 
from 101° to ~ 124° while the latitude systematically 
changed from 7^ —22° to ~ —14°. The deflection rate is pre- 
sented in Figure fToT dV It is found that the event has higher 
deflection rate at the lower altitudes, and the deflected rate 
quickly decreased to about l°/i?s beyond ~ 6 i?s- 

The deflections and the magnetic field energy density dis- 
tributions of all the data points of CME-3 are compared as 
shown in Figure [TT] The synoptic chart of the 2063 Carring- 
ton rotation, which begins at 10:03 UT, 4 November 2007 
and ends at 17:29 UT, 1 December 2007 is used. From this 
figure, it is found that the gradients of the magnetic energy 
density are well aligned to the deflections. The previous two 
CMEs occurred far away from the HCS, and then deflected 
towards the HCS. However, this CME almost initially orig- 
inated near the HCS, and deflected along the HCS. 

3.4 The 3 November 2008 Event (CME-4) 

This CME first appeared in the CORl FOV at about 23:35 
UT on 2 November 2008. The first and last images of the 
CORl data were taken at 00:05 UT and 01:05 UT on 3 
November 2008, respectively, and there are seven image pairs 
during the interval. The first and last images of the C0R2 
data are taken at 03:22 UT and 06:52 UT, respectively, and 
a total of eight data points are selected. 

Same as the CME-3 event, there is no clear source region 
observation. The tilt angle and half angle are fixed to —10° 



Table 4: The fitted free parameters derived by the model 
with the tilt angle of —10° and the half angle of 11° of the 
3 November 2008 CME. 



Time (/)c 
{VT\ \deg\ 



{deg\ \deg\ 



hi/Rs 



CORl 



00:05 




7.6 


22.6 




n 9"^ 




00:15 




7.7 


21.8 


o.ou 




< 


00:25 


234.6 


6.9 


21.3 


3.61 


0.23 




00:35 


236.6 


9.0 


21.2 


3.80 


0.23 




00:45 


236.5 


8.9 


21.6 


3.93 


0.23 




00:55 


235.9 


8.5 


20.5 


4.27 


0.23 




01:05 


236.8 


9.5 


19.6 


4.57 


0.23 








C0R2 








03:22 


236.9 


10.8 


18.2 


7.51 


0.23 


1 


03:52 


236.1 


10.3 


18.4 


8.50 


0.23 




04:22 


235.5 


9.9 


18.4 


9.41 


0.23 




04:52 


236.4 


11.2 


18.1 


10.14 


0.23 




05:22 


235.5 


10.5 


18.5 


10.95 


0.23 




05:52 


235.2 


10.5 


18.6 


11.97 


0.23 




06:22 


233.6 


9.2 


19.0 


12.64 


0.23 




06:52 


234.6 


10.4 


19.0 


13.09 


0.23 



and 11°, respectively, to get the best fitting of the CME 
shapes observed by both the STA and STB spacecraft. The 
CCS flux rope can fit the CME fairly well as shown in Figure 
1121 Table |4] lists the parameters of all the 15 data points of 
the event. 

Figure [T3j a) shows the height-time plot of the CME, and 
suggests a slight acceleration during the propagation. The 
average speed of the CME is 285.8 km s~^, and the average 
acceleration is 4 m s~^. The variations of the Stony hurst 
longitude and latitude of the CME are shown in Figures 
fTsT b) and [TsT c). respectively. Neither the Stonyhurst lon- 
gitude nor the latitude shows a significant change. There 
is no more than 3° variation of the longitude. The latitude 
just changed slightly from 22.6° to 19.6° in the CORl FOV 
and did not vary in the C0R2 FOV. Moreover, its deflection 
rate is no more than 3° / Rs as shown in Figure [TSjd). Thus, 
this CME could be taken as the event without an obvious 
deflection during its propagation. 

The comparison between the deflections and the magnetic 
field energy density distributions of all the data points are 
presented in Figure [Til The synoptic chart of the 2076 Car- 
rington rotation which begins at 23:42 UT, 23 October 2008 
and ends at 07:00 UT, 20 November 2008 is adopted as the 
bottom boundary of the CSSS model. From the figure, it 
can be seen that the gradient direction of the magnetic en- 
ergy density is toward to the nearby HCS. The changes in 
the CME propagation direction almost do not align with the 
directions of the gradient. Compared to the previous three 
events, both the deflection rate and the gradient for this 
event are quite small. Thus the large deviation between the 
two directions might not be an inconsistency. 

4 Statistical Analysis of CME Deflections 

In the above analysis, we present four events which man- 
ifest different deflection properties during their propagation 
in the corona. The first two events: the 12 December 2008 
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Figure 12: The fitting example of CME-4. From (a) to (c): The original CME images, the modeled wireframe images 
which overlays on the CME images, and the relative brightness derived from the GCS model. The top and bottom panels 
present the results based on the STA and STB data, respectively. 
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Figure 14: Same as the Figure [H but for the 3 November 2008 event. 
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event (CME-1) and the 9 April 2008 event (CME-2), which 
appeared apart from the hehospheric current sheet at the 
early stage, deflected basically along the gradient direction 
of the magnetic energy density. Both of them approached 
toward the HCS which is generally located at the region with 
the lowest magnetic energy density. The 16 November 2007 
event (CME-3), which initially originated near the HCS, 
manifested a deflection along the HCS that is in both lon- 
gitudinal and latitudinal directions. The 3 November 2008 
event (CME-4) did not exhibit an evident deflection, and ac- 
cordingly the gradient of magnetic energy density was also 
very small. 

Besides the above four events, another six events which 
have clear observations in the FOVs of the CORl and the 
C0R2 during the period from November 2007 to the end of 
2008 have also been studied. The main parameters of all the 
ten events are listed in Table [5] For each event we give the 
parameters at the first and the last valid times and list them 
in two rows. 

It could be found from Table [5] that there are two events, 
22 January 2008 and 3 November 2008 (CME-4), which 
did not deflect obviously. The other events all manifested 
a deflection during the propagation. Especially the events 
erupting on 16 November 2007 (CME-3) and 13 November 
2008 deflected in both longitudinal and latitudinal direc- 
tions. Similar to CME-3, the 13 November 2008 CME also 
deflected along the HCS. 

For all the events, we have a total of 118 data points of 
deflection, magnetic energy density gradient, and the corre- 



sponding information such as height, the instantaneous ve- 
locity, etc. Figure [15^ a) shows the distribution of the angle 
between the directions of the deflection and the gradient 
of magnetic energy density. It can be read from the his- 
togram that the fraction of events decreases from small an- 
gles (agreement between the direction of the gradient of mag- 
netic energy density and the deflection) to large angles (dis- 
agreement). The bin of the angle < 15° has the most data 
points, almost half the data points have the angle < 45°, and 
as much as 80% data points have the angle < 90° . The angle 
< 90° means that the deflections are marginally consistent 
with the gradient directions, and the angle < 45° indicates 
a good consistency. Figure [TSjb) presents the probability 
of the angle < 90° (diamond) and < 45° (asterisk), respec- 
tively, as a function of the height of the CME leading edge. 
At any height, there are at least 80% of data points with 
the angle < 90°, and at least 45% of data points with the 
angle < 45°. Particularly, the probabilities are higher within 
about 6 Rs , which suggests that the deflections and the gra- 
dients have a better consistency in the inner corona. Figure 
I15r cl presents the probability versus the strength of the mag- 
netic energy density gradient. It is clearly shown that the 
deflections and the gradients have a better consistency when 
the gradient is stronger. When the gradient is larger than 
5 X 10~^^J km~^, the two directions are marginally consis- 
tent, while the gradient is larger than 5 x 10~^ J km~^, they 
are highly consistent. 

A further quantitative analysis about the deflection rate 
is shown in the Figure 1161 The data points with the angle 



Table 5: The fitted free parameters of all the CME events. 



Date 


Time 

[UT] 


0c 

[deg] 


0s 

[deg] 




[deg] 


hf/Rs 




7 
[deg] 


a 
[deg] 


Aa 
[deg] 


V 

[km s~'^] 


N 


16-NOV-07 


09:35 
15:22 


303.1 
322.4 


101.1 
123.6 


-22.2 
-13.6 


3.0 
12.0 


0.25 
0.27 


-24.6 


8.4 


24.0 


255 


11 


04-Dec-07 


06:05 
16:52 


34.3 
34.2 


67.6 
73.4 


14.2 
3.9 


3.3 
13.8 


0.14 
0.21 


-56.5 


10.1 


11.9 


182 


19 


22-Jan-08 


22:45 
06:22^ 


226.7 
223.9 


194.6 
196.0 


-23.8 
-25.2 


3.7 
16.4 


0.27 
0.36 


-21.8 


10.1 


1.9 


296 


16 


23-Feb-08 


17:05 
06:22^ 


224.1 
211.2 


250.2 
244.6 


21.5 
19.2 


3.5 
15.5 


0.21 
0.24 


25.7 


11.7 


6.0 


174 


25 


25-Mar-08 


19:05 
20:52 


194.2 
201.7 


270.0 
278.4 


-12.8 
-12.3 


2.5 
12.7 


0.25 
0.30 


35.2 


12.6 


8.5 


1092 


4 


05-Apr-08 


16:15 
18:22 


260.8 
251.9 


120.1 
112.3 


1.3 
3.8 


3.3 
14.0 


0.17 
0.25 


-64.8 


9.8 


8.1 


982 


5 


09-Apr-08 


10:45 
14:52 


187.6 
201.6 


96.6 
112.9 


-21.9 
-18.5 


2.3 
12.7 


0.22 
0.22 


8.4 


10.6 


16.6 


476 


9 


03-NOV-08 


00:05 
06:52 


235.5 
234.6 


7.6 
10.4 


22.6 
19.0 


3.2 
13.1 


0.23 
0.23 


-10.1 


11.2 


4.5 


286 


15 


13-NOV-08 


13:05 
21:22 


288.8 
275.5 


199.9 
191.1 


-22.3 
-12.1 


3.0 
16.4 


0.24 
0.27 


-30.2 


11.2 


13.4 


256 


24 


12-Dec-08 


05:35 
14:52 


72.8 
74.3 


2.0 
8.7 


30.7 
9.6 


2.5 
17.4 


0.22 
0.29 


-15.1 


14.0 


22.1 


276 


24 



* For each event, the parameters at the first and last valid times are given in two rows. Column "Date" gives the date 
when the CME occurred. The second column lists the time when the CME observed and the superscript 'n' means the 
time of the next day. The next seven columns give the model parameters: the Carrington longitude ^(j)c\ the Stony hurst 
longitude '(/)s', latitude height '/i/', ratio ^k,\ tilt angle '7', and half angle 'a'. The 10th column means the solid angle 
between the first and the last data point. The next two columns give the average speed of the event and the total number 
of data points. 
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Figure 15: Statistical analysis of the angle between the deflection direction and the magnetic energy density gradient, (a) 
The distribution of the angle, (b) The probabilities of the angle < 90° (diamond) and < 45° (asterisk) as a function of 
height, (c) Same as the figure (b), but presents the probabilities as a function of the strength of the magnetic energy 
density gradient. 
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Figure 16: Quantitative analysis about the deflection rate. The red diamonds mark the 'bad' points, at which the deflection 
direction is opposite to the gradient direction. Three panels show the scatter plots between the deflection rate versus (a) 
height, (b) the strength of the gradient, and (c) the instantaneous radial speed. In panel (b) the instantaneous radial speed 
of a CME is coded in gray scale; darker colors stand for larger speeds. 
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Figure 17: The fitting example of CME-1 due to the GCS model with the tilt angles of —15° 
(c): The original CME images, the modeled wireframe images with the tilt angle of —15° and — i 
panels present the results based on the STA and STB data, respectively. 



and -84°. From (a) to 
)4°. The top and bottom 



between the directions of the deflection and the magnetic 
energy density gradient larger than 90° are defined as 'bad 
points' and indicated by the red diamonds. Figure llGf a) 
shows the deflection rate as a function of the height. It is 
found that the deflection rate of the CMEs decreases quickly 
with increasing height. The CMEs generally have large de- 
flection rate within about 4 i?s, and in the outer corona, the 
deflection rate approaches to zero. None of the bad points 
is beyond the deflection rate of 3°/i?s. 

The correlation between the deflection rate and the 
strength of the magnetic energy density gradient is shown 
in Figure llOf b) . The instantaneous radial speed of a CME 
is represented in gray scale, in which darker symbols stand 
for larger speeds. It is found that all the 'bad points' (shown 
in red diamonds) appear in the area with the magnetic en- 
ergy density gradient lower than 5 x 10~^J km~^ and the 
deflection rate within 3° /Rs. In such a region, any errors in 
our calculation may become relatively significant. Thus, as 
we have stated before, these 'bad points' cannot be treated 
as an inconsistency. Without these bad points, there is an 
evident positive correlation between the deflection rate and 
the gradient strength. The correlation coefficient is about 
0.85. It suggests that a stronger gradient causes a larger de- 
flection rate. Considering that the deflection rate is not sig- 
nificant when the magnetic energy density gradient is lower 
than 10~^ J km~^, we also show the correlation between the 



deflection rate and the gradient for the data points with the 
gradient larger than 10~^J km~^, which is represented by 
the solid line in Figure llGf b) . It still shows a positive cor- 
relation, though the correlation coefficient has decreased to 
0.58. 

The relative low correlation coefficient is mainly due to sig- 
nificant scatter of the data points at the right-upper corner. 
We notice that there are three data points at strong gradi- 
ents but having a relatively small deflection rate (marked by 
'A'), and other three data points with a large deflection but 
at relatively weak gradients (marked by '□'). The two sets of 
data points are obviously against the overall correlation we 
obtained, and may imply that there should be other factors 
influencing the deflection rate of CMEs. 

A possible factor is the CME radial speed. By comparing 
the deflection rate of these data points with the CME radial 
speed (derived from the height-time plot, see Figure [3ja) 
for example), we found that the speeds of the data points 
'A' are bigger than those of the data points '□'. For all the 
other data points. Figure [T6jc) shows the deflection rate as a 
function of the CME radial speed. A weak ant i- correlation 
is found between the deflection rate and the speed. The 
data points with a higher speed generally experience a slower 
deflection, i.e., the faster a CME moves outward, the smaller 
is the deflection rate. 

Besides, the CME mass should be another impor- 
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Table 6: The fitted free parameters of 
CME derived by tiie GCS model with 
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tant factor. Mass characterizes the inertia of a CME. 
Thus, the heavier a CME is, the smaller should be 
the deflection rate. However, there are only four 
events having available mass in the CDAW CME catalog 



(http://cdaw.gsfc.nasa.gov/CME_list/). The event number is 
too small to derive a reliable result. Moreover, consider- 
ing tha t there are significant errors in the mass determi- 
nation [C olaninno and Vourlidas\ "2009*; \Luaaz etal I 12005^; 
I Vourlidas et a/., 2000; Wang a/., 2011], the effect of CME 
mass on the deflection is not analyzed in this paper. 



5 Summary and Discussion 

In this paper, the deflections of ten CMEs which occurred 
from November 2007 to the end of 2008 were studied. With 
the aid of the GCS model, eight of these CMEs are found 
to be deflected during their propagation in the corona. The 
distribution of the coronal magnetic field extrapolated from 
the SOHO/MDI magnetic synoptic charts suggests that the 
CMEs tend to deflect to the region w ith lower magnetic en- 
ergy density. It confirms the result oi \Shen et al\ [201 ll ]. 

The further quantitative analysis here reveals that the de- 
flections and the magnetic energy density gradients have a 
better consistency in the lower corona or in the region with 
a stronger gradient of the magnetic energy density. The 
comparison of the deflection rate to the CME height and 
the speed suggests that CMEs have higher deflection rates 
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Figure 18: Comparison between the kinetic evolutions of the 
CME-1 event which were derived by the GCS model with 
two different tilt angles. The asterisks and the diamonds 
present the GCS model with the tilt angles of —15° and 
—84°, respectively. Panels (a) to (d) show the height-time, 
longitude- height, latitude- height, and deflection rate- height 
curves, respectively. 



in the inner corona, generally below 4 Rs. There is a pos- 
itive correlation between the CME deflection rate and the 
strength of the magnetic energy density gradient. A stronger 
gradient may cause a larger deflection rate. Meanwhile, the 
CME speed has a negative effect on the deflection rate. A 
faster event tends to have a slower deflection. It is due to 
the gradient force of the magnetic energy density acting on 
the fast event lasted much shorter than that acting on the 
slow CME. 

The fixed tilt angle implies the hypothesis that the CME 
did not rotate during the period of interest. We realized that 
this hypothesis may not be true. But it does not affect our 
result, because, even if we adjust the value of the tilt angle, 
the key parameters we have derived, including the longitude, 
latitude, and height, do not significantly change as long as 
the GCS flux rope can fit the observed CME well. Listed 
below are two examples illustrating this issue. 

In our study, the tilt angle of the CME-1 was fitted to 
— 15°. When we adjust the tilt angle, it is found that the 
GCS flux rope with the tilt angle of —84° also fits the ob- 
servations well. Figure [T7] compares the GCS model result 
between the two tilt angles, and Table [6] listed the parame- 
ters for the comparison with Table □ 

At the tilt angle of —15° the GCS flux rope is in axial- 
view, while at the tilt angle of —84° the GCS flux rope is in 
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Figure 19: The fitting example of CME-2 witfi tfie GCS model of different tilt angles. From (a) to (c): The original CME 
images, the modeled wireframe images with the tilt angle of 8° and —39°. The top and bottom panels present the results 
based on the STA and STB data, respectively. 

This event was studied in which the 

tilt angle was chosen as —53°. By comparing our results 
with the parameters which derived by the GCS model with 
the tilt angle of —53° fo r the data point r ecorded at 12:52 
UT given in the paper of iLm et al\ |2QlQbt . it is found that 
the differences in longitude and latitude are both 2° only. 

Similar to the CME-1 event, we find that the GCS flux 
rope with the tilt angle of —39° also fits the observed shape 
of CME-2. The fitting results of the GCS model with two 
different tilt angles are presented in Figure [19] and listed 
in Table [T] Although the tilt angles are different, the key 
parameters do not have significant differences. The maximal 
differences in the height, longitude, and latitude between 
the two different tilt angles are about 0.1 Rs, 2° and 2°, 
respectively. Also the evolutions in the height, longitude, 
and latitude of the CME under the two different conditions 
are quite similar, as shown by the asterisks and diamonds in 
Figure [201 

In addition, six of our ten events w ere also listed in Table 
[T] of the paper by \ Thernisien et al ] [2009]. By comparing 
our Table [5] with their table, we find that the difference be- 
tween the longitudes is mostly within 4° and four of the six 
events are just 2°, and the difference between the latitudes 
is less than 1° except for one event, which is about 3°. It 
is noticed that the longitudinal difference for the data point 
recorded at 17:52 UT of the 5 April 2008 event is about 14°. 



Table 7: The fitted free parameters of the 9 April 2008 CME 
derived by the GCS model with the tilt angle of —39° and 
the half angle of 10°. 
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side-view. Although the tilt angles are very different, the key 
parameters do not have significant differences. The maximal 
differences in the height, longitude, and latitude between the 
two cases are about 0.2 i?s, 4° and 2°, respectively. The 
variations in the three parameters are also shown by the 
asterisks and diamonds in Figure [TSl It can be seen that the 
two symbols are almost overlapped. 
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Figure 20: Comparison between the kinetic evolutions of the 
CME-2 event which were derived by the GCS model with 
two different tilt angles. The asterisks and diamonds present 
the GCS model with the tilt angles of 8° and —39°, respec- 
tively. Panels (a) to (d) show the height-time, longitude- 
height, latitude- height, and deflection rate- height curves, re- 
spectively. 



By fitting the CME im ages with the parameters given by 
\Thernisien et al\ [2009^ . we find that the difference is caused 
by the different front edge selection. Even if we adopted the 
CME front edge selected by them, and performed the same 
analysis, it can be found that the CME would manifest the 
same deflection behavior. 

In summary, the CME deflection is mainly controlled by 
the gradient of the coronal magnetic field based on our 
statistical study. The r esults confirm tha t the theoretical 
method proposed by the \Shen et al\ [201 1^ is able to quan- 
titatively describe the CME deflections. Moreover, we be- 
lieve that the method can be developed into a promising 
model, magnetic energy density gradient (MEDG) model, 
of predicting the CME deflection in the corona, though the 
basic idea of it is very simple. In this model, the gradient 
of the magnetic energy density is treated as a major cause 
of the CME deflection. Actually, the polarity of the back- 
ground magnetic fi eld may also have some effect on the de- 
flecti ons of CMEs \Chane et all l2005l : \lsenherQ and Forhe^ . 
12007^ . Besides, it should be noted that the gradient of the 
magnetic energy density decreases quickly with increasing 
height. When a CME propagates outward, the gradient of 
the background magnetic field may become weak rapidly. 
Such weak gradient would not be sufficient to make a CME 
deflected obviously, particularly during the propagation of a 



CME in the interplanetary space. This implies that there 
should be another mechanism to cause the CME deflec- 
tio n in the interplanetary spa ce, which had been reported 
bv lPoommses et al\ [2010[ | cind \LuQaz et al\ [2010l ]. A possi- 
ble candidate mechanism is the CME's i nteraction w i th the 
background solar wind as proposed by I Wang et al. 1 [2004 
2006b]. 
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